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Q_i. Abstract 
< 

We critically re-examine the case for and against a sizeable right-handed component 
in the b — ► c charged current coupling with a strength £ relative to the conventional 

> 

in 
oo 

■ knowledge of Uft extracted from both inclusive and exclusive processes, we are able to 

m ; 

determine the presently allowed parameter space for £ via HQET. We then identify 

OO 

several observables which could be measured at B factories to either strengthen these 
Q_i| constraints or otherwise observe right-handed currents. This parameter space region is 

found to be consistent with the low degree of A& polarization as determined by ALEPH 
as well as the measurements of the charged lepton and neutrino energy spectra from 



left-handed current. Using data from CLEO on the decay B —* D*£v, as well as our 



^ ■ b decay made by L3. We discuss how future measurements of semileptonic decay 

distributions may distinguish between exotic Af, depolarization mechanisms and the 
existence of right-handed currents. Within the parameter space allowed by CLEO, 
using the Left-Right Symmetric Model as a guide, we perform a detailed search for 
specific sub-regions which can lead to a reduction in both the B semileptonic branching 
fraction as well as the the average yield of charmed quarks in B decay. The results 
provide a concrete realization of an earlier suggestion by Voloshin but may lead to 
potential difficulties with certain penguin mediated decay processes. 
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1 Introduction 



The new generation of B factories soon to come on line will open up an additional window 
on potential physics beyond the Standard Model(SM). Unlike the situation at new higher 
energy colliders, the physics beyond the SM will appear indirectly, e.g., as deviations from 
SM expectations in decay rates, distributions and/or asymmetries obtained through precision 
measurements. Such measurements may be as important in probing the SM as are those 
currently performed by LEP/SLC and the Tevatron at higher energies and will rival others 
associated with the observation of CP violation in the B system. 

Perhaps the most fundamental of all quantities associated with b quark decay is the 
chirality of its charged current coupling. The possibility that the b — > c charged current (CC) 
may have a sizeable right-handed (RH) component has been the subject of speculation for 
some time. Early on, Gronau and Wakaizumi, as well as a number of other authors 
speculated that the b — > c coupling might be almost, if not purely, RH. Thanks to mea- 
surements performed by both the CLEOp| and L3|§ Collaborations, to which we return at 
some length below, we now know that this hypothesis cannot be true. The relative strength 
of the RH b — > c coupling in comparison to the corresponding SM left-handed(LH) coupling 
must be somewhat less than unity; the leptonic current in the decay is highly constrained 
to be LH. It is important to observe that the results of these experiments cannot exclude 
a RH coupling of modest strength. Other data, such as the apparently small value of the 
Ab polarization observed in Z decay[|J by ALEPH, qualitatively support the hypothesis of 
potentially sizeable RH couplings unless some exotic depolarization mechanisms are at work. 
Thus the current experimental situation remains unsatisfying and is far from resolving the 
issue of the presence of RH couplings in b — > c transitions. On the theoretical side, in a 
completely different context, Voloshin|| has recently suggested that a RH b — > c coupling 
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of modest strength may help to resolve the well-knowncite blok B semileptonic branching 
fraction(I^) and charm counting(n c ) problems [|7|]. Thus we are left with the questions: given 
the present data, it is possible that such a RH coupling actually exists and can it assist with 
the Bi and n c problems? How can future B factory data help clarify this situation? 

In this paper we will examine the simultaneous compatibility of the CLEO and, to 
a lesser extent, the L3 constraints on the RH b — > c coupling and the desire to address 
the charm counting/branching fraction problem along the lines suggested by Voloshin. The 
important role played by the ALEPH polarization measurement is also examined. This 
discussion stresses both what we can learn from the current data about possible b — > c RH 
couplings and what can be learned through future precision measurements at B factories to 
resolve the present ambiguous situation. Most of this analysis, associated with the constraints 
from the CLEO, L3 and ALEPH data, can be performed in a completely model-independent 
fashion without any reference as to the possible origin of the b — > c RH coupling. However, 
in order to subsequently approach the n c — Bi problem a more rigid theoretical framework, 
such as the Left-Right Symmetric Model (LRM)||, needs to be invoked. Other more general 
frameworks are possible but are beyond the scope of the present paper. 

This paper is organized as follows. In Section 2, we re-examine and update the 
constraints imposed on the relative strength of the RH to LH b — > c coupling due to current 
experimental data from CLEO, ALEPH and L3 using a model-independent approach and 
relying on Heavy Quark Effective Theory. As we will see there is a tantalizing, though 
still not compelling, hint of RH interactions in the CLEO data. The importance of both 
the small observed A& polarization as a possible, though still ambiguous, signature for RH 
currents is discussed. We also examine how this scenario may be distinguished from the SM 
with exotic depolarization mechanisms. The present ALEPH data is shown to be consistent 
with a moderately strong RH current coupling. In examining how we can extract further 
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information from the present data and with an eye towards future measurements we discuss 
several new observables and their potential usefulness in probing for RH couplings. Many 
of these observables have either not been measured or have yet to be examined with any 
degree of precision. Several of these quantities can be probed at the Z or during the first 
year of running at the new B factories. In Section 3, we present an overview of the LRM 
and discuss the meaning of the current experimental constraints within this specific context 
keeping detailed discussion of the required LRM particle content to a minimum. As we do 
not want to restrict or constrain ourselves to a specific version of this model, we tacitly avoid 
at this point any discussion of loop processes which may involve the full particle spectrum of 
a realistic, probably supersymmetric LRM. In Section 4 we describe the nonleptonic b — > c 
decays in the presence of RH currents and their associated decay widths including LO and 
estimates of the NLO QCD corrections based upon what is currently known in the case of the 
SM. In Section 5, we use the LRM as input to scan the full model parameter space allowed 
by the CLEO data to discover and identify sub-regions that will lead to a simultaneously 
decrease in the values of both and n c in comparison to the SM expectations. While 
such regions are shown to exist they occur in only a small, fine-tuned, fraction of the entire 
parameter space volume. Our summary and conclusions can be found in Section 6. In 
the Appendix we speculate on the possible forms of Vr and point out how certain penguin 
processes can lead to difficulties with the solutions to the n c —B^ problem that we've obtained. 

2 Constraints on Right- Handed b — > c Couplings 
2.1 Model-Independent Notation 

Allowing for both LH and RH b —>■ c couplings while, following Voloshin, maintaining the 
leptonic current as purely LH to satisfy the well-known /i decay constraints^] without fine- 
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tuning neutrino masses, the general four-fermion interaction describing B semileptonic decay 
can be written as 



AC 

7F -V^ b [{c Ll ,b L )+i{c Rl ,b R )]{I Ll ^L) 



v. 



V2 

where here we will treat £ as a complex parameter, £ = |£|e* A , though CP violation will 
be ignored in the discussion that follows [|Kj. As we will see the additional phase degree of 
freedom will play a very important role in obtaining signals and constraints on RH currents. 
We recall that in the original Gronau and Wakaizumi scenario the leptonic current in B 
decays was also RH|J and neutrino masses were tuned to allow for the semileptonic decay 
process. How do we ascertain the allowed range of £? Again following Voloshin, the first 
place to obtain a constraint is inclusive semileptonic b decay at the quark level. The most 
obvious observable is the inclusive decay partial width that can be written as 



r(b^c£is)~\v£\ 2 f(x)ri 



le| 2 + 2i?e(0 



9{x) Vr 



(2) 



For zero mass leptons, f,g are the well-known kinematic phase space functions]!, |TT], [12| of 
the ratio x = m c /mb ~ 0.29: 



/ = (l-x 4 )(l-8a; 2 + a; 4 )-24x 4 lnx, 

g = -2x[(l -x 2 )(l + 10x 2 + x 4 ) + 12x 2 (l +x 2 )\nx] 



(3) 



and for x = 0.29 we find / ~ 0.542 and g ~ —0.196. For semileptonic decay to r's, the 
corresponding phase space suppression factors can be decomposed as f T = {I\ + J 2 )/2 and 
g T = (I i — -?2)/2 in terms of the integrals 



h = ds Z 

Jv 2 



ds Z 



A(4s - y 2 ) + 2AE(1 + 2y 2 /s) - (E + 2s)(2s + y 2 



E(4s - y 2 ) + 2AE(1 + 2y 2 /s) - (A + 2s)(2s + y 2 ) 



(4) 



where y = m T /mb, £ = (1 + x) 2 and A = (1 — x) 2 with 



y 2 

1 - — 



2 r 



(S-a)(A-a) 



1/2 



(5) 



Numerically, one finds f T ~ 0.122 and g T ~ —0.0490 for x = 0.29 and ?/ ~ 0.372. Note that 
when RH currents are present the ratio of the b semileptonic decay width to r's to that for 
massless leptons becomes a weak function of £ with overall variations of order a few per cent; 
this dependence occurs due to a mismatch in the phase space ratios: g T / f T 7^ g/f- This 
effect is most likely too small to be observed experimentally, however, but should be kept in 
mind. 

The parameters tj^r represent both perturbative and non-perturbative strong inter- 
action corrections which also depend on x as well as the lepton mass and the relevant strong 
interaction scale //. When needed in the numerical discussion below we will assume that 
the effects of all strong interaction corrections in the b — > c semileptonic decay are at least 
approximately insensitive to the chirality of the charged current coupling, i.e., T]l = Vr = V 
as was done Voloshin. (Certainly, an explicit calculation needs to be performed to verify 
this assumption.) To leading order in QCD for massless leptons and with x = 0.29 one 
has the perturbative contributions r\ = 1 — ^(2.53) + 0(a 2 ) whereas for r's one obtains 



rj — 1 — ^(2.11) + 0(a^)|13fl. The complete NLO expressions are not yet available in either 



case. Only the terms of order a 2 j3 , with (3q being the one-loop QCD beta function, are 



known at present |I1 ] , so for now we will truncate these corrections at this order but include 
them in our detailed numerical analysis below. 

It is amusing to note that the existence of a RH coupling means that a measurement 
of the partial width Y yields not the true but an effective value of from inclusive data 
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when the result is interpreted in terms of the SM, i.e., 



\TrL\inc IT/^I 

I V cb\eff — \ V cb\ 



1/2 

(6) 



under the assumption that t)l = Vr- This result will have important consequences for us 
below. 

To obtain more information from this inclusive decay additional observables are re- 
quired. Indeed, many authors have speculated about how one can experimentally extract 
information about potential RH couplings in inclusive semileptonic b decay. Dittmar and 



Was|15[ suggested examining simultaneously both the charged lepton and neutrino, i.e., 
missing energy, spectra arising from b semileptonic decay at the Z peak. When one looks 
at the squared matrix element for this process in the free-quark limit, after all traces are 
performed, the sensitivity to £ becomes immediately transparent: 

\M\ 2 ~ pt ■ p c p u ■ p b + \£\ 2 pe ■ p b p u ■ p c - m h m c pi ■ p v Re(£) , (7) 

with the pi labelling the corresponding particle four-momentum. The £ sensitivity is seen to 
be particularly enhanced due to the large value of the ratio m c J m& ~ 0.29 with the phase of £ 
playing a very important role. (For completeness we note that one can find the full expression 
for the resulting unpolarized charged lepton spectra in the Z rest frame at leading order is 
given by Fujikawa and Kawamoto ||1 1||. The corresponding neutrino spectrum can be trivially 
obtained through the interchange of the LH and RH couplings.) Interestingly, as mentioned 
earlier, L3|| performed a simultaneous measurement of both the charged lepton and missing 
energy spectra in b decay and excluded very large values of £, i.e., purely RH couplings, by 
more than 6cr and £ ~ 1, i.e., purely vector couplings, by more than 3cr. They did not, 
however, attempt a fit to £ as the required sensitivity to £ << 1 was not available once 
detector cuts and hadronic as well as other systematic uncertainties were taken into account. 
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However, values of |£| << 1 were certainly not excluded and we will attempt to further 
quantify these results below. 

We now turn to each of the three experiments CLEO, ALEPH and L3 and survey the 
constraints that they are presently imposing on £ and what can be learned from comparable 
measurements at future B factories even if relatively low integrated luminosities are available. 



2.2 CLEO 

In addition to inclusive semileptonic decay one may hope to obtain information on possible 
RH coupling through exclusive decay measurements due to the enriched nature of the acces- 
sible final states. In this regard CLEOjpJ has performed a detailed examination of both the 
B — > D and B — > D* exclusive semileptonic modes. In the B — » D case the impact of RH 
currents is well known to be rather minimal for massless leptons since the final state and 
the corresponding hadronic matrix element are rather simple. In this care, their only effect 
is to scale the anticipated partial rate by an overall factor, |1 + £| 2 , to which we will return 
below. A more complex and interesting pattern occurs for the B — > D* case. 

The CLEO analysis^] that examined the exclusive decay B° — > D~*(^ Dtt)£u sought 
to extract form factor information and, in particular, to measure the forward-backward 
asymmetry of the charged lepton, Afb, the average D* polarization, Tl/Tt, as well as 
Vj£. The data sample of ~ 780 events employed in their analysis resulted from an initial 
set of 2.6 ■ 10 6 BB's corresponding to an integrated luminosity of ~ 2 Afb" 1 at the T(AS). 
Following the general analysis as presented in Ref . 0, |16| , one begins with an initially four-fold 
differential distribution but this is a bit unwieldy. Integration over two of the three decay 
angles (the others of which we will subsequently return to below) leads to the following 
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double differential decay distribution for this process in the massless lepton limit: 



d 2 Y 

dq 2 dz 



~ \V^\ 2 Pq 2 



(1 - z) 2 \H + \ 2 + (1 + zf\H4 2 + 2(1 - z 2 )\H \ 



(8) 



where P is the D* momentum in the B frame, q 2 is the four-momentum transfer from the B 
to the D* and z = cos 0g with 9g being the decay angle of the i in the virtual W rest frame. 
P is given by 



P 



2M 



(M 2 - m 2 - q 2 ) 2 - Am 2 q 2 



1/2 



(9) 



and the helicity amplitudes H±$ are functions of q 2 which are generally expressed in terms 
of the conventional form factors A\ 2 and V as 



H ± {q 2 ) = (M + m)A 1 (q 2 ) + 



2MP 
(M + m) 



V{q 2 ) 



H (q 2 ) = [2m J q 



21-1 



4 /If 2 P 2 

{M 2 -m 2 - q 2 )(M + m)A l (q 2 ) - A 2 (q 2 ) 

(M + m) 



(10) 



where M{m) is the mass of the B(D*). Meticulously following Neubert|T7[ one may use 
suggestive versions of the above form factors that have very well defined limits when Heavy 
Quark Effective Theory (HQET) becomes exact: 



Mq 2 ) 

V(q 2 ) 



(M + m) 
2\[Mm 

(M + m) 
2\[Mm 

(M + m) 
2VM> 



m 



(M + m) 2 
R 2 (w)h(w) , 



h(w) 



(11) 



Here we define as usual w = (M 2 + m 2 — q 2 ) / (2Mm) . In the exact HQET limit both R12 — > 1 
and h(w) becomes the Isgur-Wise function so that the Ri can be considered as representing 



small corrections in both a s and 1/m to the case of pure leading order HQET. Generically, 
h has a linear form, h(w) = h(l)[l — p 2 (w — 1)] although other structures are possible. 

While the forward-backward asymmetry can be obtained by integration of the expres- 
sions above, the ratio Tl/Tt can be determined from the decay angular distribution of the 
D in the D* frame when D* — > Dti (cos6V, in the notation of Ref.0). Following Ref.[0, [16 
we can write the relevant double-differential distribution in this case as 



d 2 T 



dq 2 d cos 9 V 



Pq 2 



(\H + \ 2 + |#-| 2 )(1 - cos 2 6 V ) + 2\H \ 2 cos 2 6 X 



(12) 



T L /T T essentially probes the relative weights of the H and H± helicity amplitudes as we 
will see shortly. 

So far this discussion has been quite general. To include the effects of £ 7^ in 
comparison to SM expectations we simply make the replacements V — > V(l + £) and — > 
^1,2(1 — in the expressions for the helicity amplitudes above and recall that £ is complex. 
This follows directly from the rescaling of the LH and RH current amplitudes as seen in 
Eq.(l). Once particular expressions for R\ y2 and h are assumed we may directly calculate 
Afb, Fl/Tt, as well as the total decay rate, which then gives us V c ^ exc (D*). We obtain 



.4 



FB 



fdq 2 fl°dz 



d 2 r 

20 ""^ dq 2 dz 



Jdq 2 2z (l-z 2 /3)Pq 2 H 2 
Jdq 2 z (l + z 2 /3)Pq 2 (H 2 + + H_ 



2\ ' 



(13) 



where z (q 2 ) expresses a potential minimum lepton momentum cut used to identify the event: 



z = mm 



4Mpi — M 2 — q 2 — m 
2PM 



(14) 
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CLEO, for example, employs a typical lepton momentum cut of ~ 1 GeV. These expressions 
can be re-written to clearly display their £ dependence as 



A 



FB 



(1 + \^)A - 2BRe(0 ' 



4 [1 + |£| 2 - 2Re(Q]D 
3{l + \£\ 2 )E + 2FRe{0 ' 



(15) 



Experimentally, CLEO[|] finds A FB = 0.197 ± 0.037 and T L /T T = 1.55 ± 0.29, which are 
of course both consistent with SM/HQET expectations. Here A — F are a simple set of 
numbers which result from performing the double integration over the relevant kinematics. 
For a fixed set of and h, the values of A — F are completely determined subject to 
experimental cuts, and these results can be combined to constrain £. In addition, from the 
expression for the overall partial width we also obtain 



when the value is again interpreted in the SM; this result should then be compared with 
Eq.(6). Note that since one finds that —B/A ^ gj f, the apparent values of extracted from 
exclusive B — > D* and inclusive measurements will be different when Re(£) ^ 0. Demanding 
that the true take on the same value in both cases imposes an extra constraint on £. In 
order to employ this additional constraint we use the specific numerical results as provided 
in the recent review of both inclusive and exclusive semileptonic decay data by Buras|l8| to 
obtain Vj^ exc {D*) /V^ inc = 0.967 ± 0.105. This value is completely consistent with unity, 
as anticipated, but will still provides an additional requirement on £. A similar situation, 
as mentioned above, occurs in the case of B — ► D semileptonic decays where we now would 



B 
A 



nl/2 



V%\TMD*) = \V%\ ■ 1 + |£| 2 - 2ite(£) 



(16) 
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find simply 



M 2 

rL I exc , 



I W/P) = \v£\ 



l + |^ + 2i?e(£) 



17) 



which should be compared with that from the D* mode above. Given the present exper- 
imental situation|| adding this additional constraint will not influence the results of the 
fit obtained below. However, future measurements may make this an important input into 
analyses of RH currents. 

Our procedure is the following: for a fixed set of Ri^iw) and h(w) we calculate the 
integrals A — F and then perform a simultaneous x 2 fit to the CLEO results on Apg and 
T L /T T as well as to the ratio exc (D*) jV^ mc treating |£| and ca = cos A as free parameters 
(recall, A is the phase of £). Possible correlations are ignored. We then choose another set 
of i?i 5 2 and h and repeat the process. Each repetition will thus generate a 95% CL allowed 
region in the ca — |£| plane. To be specific we employ forms of Rx^{w) and h{w) suggested 
by Neubert|17j and by Close and Wambach(CW) flT9|| as well as several other sets suggested 



by the first paper in Ref. 0. As a typical example, with R^ w = 1.15[1 — 0.06(w — 1)], 
RCW = .91[1 + 0.04(w - 1)] and p 2 = 0.91 we obtain A ~ 0.116, B ~ 0.105, C ~ 0.024, 
D ~ 0.060, E ~ 0.056, and F ~ 0.044. These values do indeed reproduce the well known 



SM expectations[17] in the £ — > limit. 



The results of this fit are shown in Fig. |T| which displays the 95% CL upper bound 
on |£| as a function of ca for several different choices of Ri^ and h. The most important 
features of these results to notice are: (i) the bounds we obtain are not very sensitive to the 
exact choice of these HQET functions and (ii) the constraints on |£| are strongest when £ is 
real. We note that Voloshin's preferred range of values of £ = 0.14 ± 0.18 lie mostly inside 
the allowed region. It is clear that at the moment the existing constraints on £ are quite 
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Figure 1: 95% CL allowed region (below the curves) in the |£| — ca plane obtained from a 
fit to CLEO data as well as the experimental value of the ratio V c f exc (D*)/V^ mc . Each of 
the six curves corresponds to a unique choice of R lj2 and h. The SM limit lies along the 
horizontal axis at |£| = 0. The locations of the six x 2 minima are also shown for completeness 
and are seen to be reasonably clustered. 
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poor and that values of |£| of order 0.25 are certainly allowed by current data. We note 
that for the six sets of HQET functions used in this analysis the resulting best fit values 
for £ are reasonably clustered and indicate a magnitude ~ 0.20 — 0.35 and a sizeable phase. 
With the far larger data sets soon to be available from B factories it is quite important 
for this analysis to be be revisited and refined in the not too distant future. We note that 
a somewhat smaller allowed region results if the unpublished results from CLEO that now 
include the charged B decay modes are usedpHI. 

One might ask if there are other observables associated with this exclusive decay 
that could allow for some additional sensitivity to RH current interactions [|16j. To this end 
we briefly examine both the q 2 and \ distributions which can be measured using the re- 
coil momentum of the D* and identifying the angle between the W and D* event planes, 
respectively. Once integrated over all other variables, the deviations in both these distribu- 
tions from the SM expectations are found to be totally controlled by the value of the ratio 
A = 2|f |c A /(l + |£| 2 )- The form of the q 2 distribution can be obtained immediately from the 
double differential expression above. Fig. |2|, where we have used the HQET functions, 
of Neubert ||17|| and those of Close and Wambach||19||, shows that the normalized distribution 



is only very weakly dependent on the existence of RH currents. Specifically, we see a direct 
comparison of the SM distribution, A = 0, with that expected for the cases of A = ±0.5. 
From the figure it appears unlikely that the shape of the q 2 distribution will yield any useful 
information on RH currents unless very high precision data is obtainable. Note there is little 
difference between the curves generated with the two different sets of HQET functions. 

In the case of the normalized x distribution, the shape is controlled by a single 
parameter if all the other variables have been completely integrated over, i.e., 

^ = i(l-ncos2 X ), (18) 
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0.1 0.2 0.3 0.4 

X 



Figure 2: Normalized q 2 distributions for the process B — > D*iv. Here x = q 2 /M 2 and the 
curves correspond to the SM(solid) and A = 0.5(— 0.5) (dotted and dashed, respectively). A 
possible p t cut on the charged lepton momenta has been ignored. Results are shown for both 
the Neubert as well as the Close and Wambach HQET functions corresponding to the pair 
of curves for each case. 
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where 



_ Jdq 2 Pq 2 Re(H* + H_) 

Jdq 2 Pq\H 2 + + H 2 _ + H%)' { } 

which we may rewrite to show the A dependence explicitly as 

n = (T2 + TlA) (20) 

(2r 1 + r3) + (2T 2 -T3)A' lzUj 

with the Tj being a set of kinematic integrals. In the SM one finds that Q ~ 0.175(0.192) 
using Neubert(CW) HQET functions. Note that if instead only the even(odd) values of 
cos 9y are integrated over, the normalized x distribution picks up an additional term of the 
form 

dN dN 3^ 
dx dx 8 

where 

_ Jdq 2 Pq*ReH*(H + -H_) 

fdqi Pq%Hl + H 2 +HZ) ■ ^ 

For this variable the £ and ca dependencies become are somewhat more complex and cannot 
be expressed simply through the parameter A; £ is expressible as 



Ta 



(2T 1 + r 3 ) + (2T 2 -T3)A' 

with T4 being another kinematic integral. In the SM one finds that S ~ — 0.25(— 0.22) for 
Neubert(CW) HQET functions. 

Fig.|3| shows that Q is quite sensitive to positive values of A so that one may hope to 
get a reasonable sensitivity to RH interactions if Q could be precisely measured. Present 
data from CLEO is found to be consistent with the expectations of the SM for Q but the 

15 



Figure 3: Value of the f2 parameter which controls the shape of the x distribution as a 
function of A. The results are shown for both the Neubert (solid) as well as the Close and 
Wambach(dashed) HQET functions which are seen to yield quite similar results. 
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statistics are still rather poor. To get an idea of the potential sensitivity we have performed 
a straightforward two parameter (normalization and Q) fit to the existing binned data as 
presented in Ref.0. To obtained improved statistics in this first fit we have combined 
the data in both the cosOy > and cos 9y < regions. Unfortunately, the resulting the 
distribution of the data shows little sensitivity to Q. After background subtraction this 
fit yields Q = 0.126 ± 0.120 at 95% CL which is certainly consistent with the SM. This 
constraint subsequently implies that A lies in the 95% CL range — 3.3(— 2.8) < A < 0.71(0.75) 
for Neubert(CW) HQET functions using the results in Fig.^. As one would expect from the 
low sensitivity to negative A, our bound in this case is rather poor. 

A second more hopeful possibility is to fit the shape of the x distribution for both Q 
and £ by treating the cos 9y > and cos By < regions independently; here there is a loss of 
statistics but a dramatic increase in sensitivity to RH couplings. Following the same analysis 
as above we arrive at the results presented in Fig.|] for the allowed region in the ca — |£| 
plane for Neubert, CW as well as Isgur and Wise[^TJ HQET functions. Note that the allowed 
region resulting from this fit is somewhat sensitive to the HQET Ri choice, quite unlike the 
other observables that we have examined up to this point. Although this result seems to 
support the possibility that RH currents may indeed be present one must be hesitant to form 
such a hasty conclusion without further analysis. First, the only believable fit of this kind 
must be performed by the CLEO Collaboration and we note the apparent strong sensitivity 
of our result to the choice of the Ri HQET functions. However, it is certainly most clear 
that our understanding of potential RH currents in b decay would very much profit from 
higher precision measurements of the x distribution. This seems possible during the first 
year of T(45) running of BABAR and BELLE since the CLEO data sample used in this 
analysis corresponded to only 2.6 million BB pairs. We note in passing that using the still 
unpublished CLEO data from the charged B decay mode|20| already strengthens the case 
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for right-handed couplings based on the fit to the x distribution. 

Another question one might ask is what the allowed range for the parameters fl and 
E are given the CLEO constraints we have extracted from the earlier fit. To obtain such 
results we need to scan the |£| — ca region below the envelope of curves shown in Fig.l to get 
the extrema. We find 0.053 < tt < 0.207 and -0.345 < £ < —0.115 for the Neubert HQET 
functions; correspondingly, for the CW HQET functions we obtain 0.089 < f2 < 0.218 and 
-0.310 < E < -0.106. 




Figure 4: 95% CL fit to the shape of the CLEO x distribution assuming Neubert (dotted), 
CW(dashed) or ISGW(solid) HQET functions. The allowed region is either below the dotted 
line or within the dashed or solid enclosure. As before the diamonds locate the x 2 minima 
for the three sets of Ri. 
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As a final note, if the r polarization in the decay B — > D*tv can be measured, 



Wakaizumi has shown ||22|| that it provides yet another quantity which is fairly sensitive to 
£ 7^ 0. This decay mode will thus yield even more observables which can be used to probe 
for b — > c RH currents due to the addition of finite mass terms associated with the r. Of 
course for this mode there is a loss in statistics due to the additional phase space suppression 
due to the r mass as well as the associated r reconstruction efficiency to be dealt with. An 



analysis of these prospects is, however, beyond the scope of the present workJlO] 



2.3 ALEPH 

Unfortunately, other data cannot at present improve significantly upon the CLEO bounds 
without further employing some rather strong assumptions. For example, in principle the low 
A;, polarization observed in Z decay Q by ALEPH can be used to obtain such a constraint. 
We recall that a b quark produced at the Z in the SM is highly polarized, i.e., P = —0.935 
and radiative effects have been shown to reduce this value only slight ly|23|]. During the 



hadronization process some of the memory of the original b polarization is lost but it had 
been anticipated that in the b — ► A& process a large part of the original polarization would 
be kept |[24|| . Falk and Peskin[^4| estimated on the basis of HQET that the resulting At 



polarization would be P = —(0.69 ± 0.06). 

The ALEPH analysis is based on ~ 3 • 10 6 hadronic Z decays which yielded a sample 
of 462 ± 31 A& candidates. The method used by ALEPH to extract the value of P for the 
A;, was first suggested by Bonvicini and Randall [p5fl who noted that the ratio of the average 
values of the neutrino and lepton energies in semileptonic B decay, y =< Eg > / < E v >, was 
particularly sensitive to the polarization of the b quark. This variable, being an energy ratio, 
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is quite insensitive to b fragmentation, detector acceptance and reconstruction effects as well 
as the uncertainties in the ratio m c /mb. We note that the direct comparison of the average 
charged lepton and neutrino energies from b quark decays with theoretical expectations, as 
was done by L3||, was found to lead to substantial fragmentation uncertainties although 
values of £ of order unity were clearly excluded; we will return to the L3 data below. It 
has also been found that a s and 1/ml corrections |27| to the parton level expectations for 
y are quite small (and, hence, will subsequently be ignored). Of course these results have 
only been explicitly demonstrated in the case of purely LH couplings. In our analysis we will 
make the reasonable assumption that they remain true when both LH and RH couplings are 
present. 

The averages of Eg and E v can be calculated directly from the decay at rest spectra 
through the boost relations Ei = j(E% + /3p* L ), etc., with /3 ~ 1 and p* L being the lepton's 
momentum in the boost direction. In order to remove selection cut and energy reconstruction 
errors which produce a bias in y, ALEPH instead determined the ratio of ratios R y = 
ydata/yMc(fy where 2/mc(0) is the y value obtained from a Monte Carlo simulation employing 
the SM in the limit of zero polarization. The value of R y was then compared with the Monte 
Carlo-corrected SM theory prediction to extract the value of P. What ALEPH found was 
R y = 1.10±0.13 (including systematic errors in quadrature)that then yielded the intriguingly 
small value P = — 0. 23^23, which is significantly smaller in magnitude, by ~ 2cr, than the 
expectations of Falk and Peskin. 

To investigate the double ratio Ry in the case when RH currents are present, we must 
return to the normalized double-differential charged lepton decay distribution. In the b rest 
frame to leading order and neglecting the lepton mass we find 

dN 
dzd cos 9 
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R{x } z) + P cos 6 Q(x } z) 
[l + mf{x) + 2Re(i)g{x) 



(24) 



where z = 2E^/m& and 9 is the angle between the b and i momenta with f(x) and g(x) given 
above. Explicitly, we find that R = R LL + Rrr\£,\ 2 + ^Re{i)RRL and Q = Q LL + Qrr\£,\ 2 + 
2Re{Z)Q RL with 
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and 



Qll 
Qrr 
Qlr 



z 2 (l — x — z) 



6z 2 (l-x 2 -z) 2 
(1^1 



;i -z){l -2z + x 2 ) -2x 2 



6xz 2 (l — x 2 — z^ 2 
2(1 -z) 



(26) 



These results confirm those obtained by Tsai|28[ long ago in a different form and context. 
The corresponding expressions for the neutrino spectrum can be obtained from the explicit 
relations above by interchanging the role of the left- and right-handed labels. Using these 



results we can calculate y following Bonvicini and Randall |[25||, rescale this value by the 
SM result assuming P = 0, and include the Monte Carlo corrections of ALEPH. Given an 
assumed value for P we can then fit to the ALEPH data to obtain an allowed region in the 
|C| ~~ c a plane. The result of this analysis is shown in Fig. [| and is compared to the CLEO 
allowed region obtained above assuming the estimate of the polarization retention of Falk 
and Peskin, P = — (0.69±0.06), is correct. Here we see that at the 95% CL almost the entire 
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plane is allowed except for a possible small region on the lower right which only appears in 
the case of P = —0.75. As P increases in magnitude, we note that the allowed parameter 
space region shrinks somewhat in size. We also see from the figure that the location of the 
best fit is quite sensitive to the assumed true value of the polarization. We note that if 
there are additional dynamical mechanisms^] which could lead to a further reduction in 
the expected value of P in the SM and they could be reliably trusted quantitatively, then 
the limits we would obtain on RH b — > c couplings might be improved. In the future if the 
central value obtained by ALEPH was verified and the errors were reduced by a factor of 
two the size of the allowed region would shrink substantially and form a band approximately 
5\£\ ~ ±0.25 wide on either side of the best fit points shown in the Figure. 

If the apparent low value of the polarization as measured by ALEPH is verified by 
future experiments then there are only two conclusions. Within the SM framework there must 
be some new source of depolarization and indeed P ~ —0.23. Alternatively, right-handed 
currents are present and the true value of P is closer to the HQET expectations of P ~ —0.69 
but appears low when interpreted in terms of the SM. As discussed above, a reduction in the 
ALEPH error by a factor of two, assuming the same central value, would clearly define a small 
allowed region when combined with the results from CLEO. Unfortunately, a measurement 
of Ry alone, no matter how precise, will be able to eliminate the possibility of some exotic 
depolarization mechanism and allow us to conclude that RH couplings exist. However, an 
analysis of the higher y moments or other possible distributions may be most helpful as 
suggested by Diaconu et al. |E3] in an important paper. For simplicity we first consider only 



the ratio of the second moments of the decay distributions here. (We have not examined 
moments higher than second.) Within the SM if x = 0.29 and P = —0.23 we can uniquely 
predict the value of the quantity R2 y = 2/2/2/2(0) = 1.181, where 2/2 =< E\ > / < E% >, 
although a s and l/mf corrections are somewhat larger here. In the case where RH currents 
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Figure 5: Comparison of the envelope of the 95% CL allowed regions obtained with CLEO 
data (solid curve) with those obtainable from the ALEPH A b polarization results assuming 
P = —(0.69 ±0.06) corresponding to the dotted, dashed and dash dotted curves. The region 
below the slightly tilted horizontal curves and outside the 'nose' on the lower right-hand side 
for the case of P = —0.75 are allowed. The location of the corresponding x 2 minima for 
P=— 0.63, —0.69, and —0.75, respectively, are also displayed from right to left. 
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are present and P = —0.69, we can invert the R y relation and find |£| as a function of ca and 
then calculate the corresponding value of R 2 y We find in this case that for the central value 
R y = 1.10 one obtains 1.198 < R 2y < 1.227 apart from the above corrections; note that this 
range does not overlap with the SM expectation but the separation between the two is quite 
small. It would thus appear that simultaneous very high precision measurements of both R y 
and R 2y , as well as possible higher moments, are required in order to be able to resolve the 
ambiguity and determine if RH currents are indeed present in semileptonic b decays. Given 
the current and anticipated sizes of the errors a determination of at least these first two 
moments alone will not necessarily prove useful. 

As a second possibility we note that Diaconu et al. also suggest a number of other 
variables which can be used to probe the A b polarization. One of these is the difference in the 
charged lepton and neutrino rapidities, Ai] — rjt — i] u , where these rapidities are measured 
with respect to the boost direction. This quantity is directly proportional to the polarization 
and, being a rapidity difference, is fortunately insensitive to fragmentation uncertainties. We 
find that with RH currents contributing one obtains 



so that in the SM for P = -0.69(-0.23) we would obtain Ar] = 0.436(0.145). In the case 
of RH currents, repeating the above procedure to find |£| as a function of ca from the data 
on R y we are led to the prediction that Arj = 0.238 — 0.257, assuming that P = —0.69, 
which is quite different from either the SM expectation with a low value of P or the HQET 




(27) 





(28) 
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SM prediction. Again it appears that the RH current and exotic depolarization mechanism 
possibilities may be separable using precision measurements. However in this case we note 
that the required level of precision for this variable is far less that that for Ri y giving us 



some hope that such a separation may indeed be possible at future B factories [26 



2.4 L3 

Following the same approach as above one might attempt to further quantify the L3|| 
constraints on £ by constructing the y values using the results presented in their Table 2 and 
including some corrections associated with their Monte Carlo. This would then be similar 
to the ALEPH analysis but now one is actually probing the initial b quark polarization 
about which there is far less uncertainty. Of course, in principle, only L3 can perform this 
procedure but our rudimentary study will provide an indication for the location and size of 
the allowed region associated with their data. If we simply double their errors but then ignore 
both the a s and I /ml corrections as well as fragmentation and energy scale uncertainties 
and neglect any correlations, we can obtain an estimate for the associated allowed region in 
the ca — |£| plane. This most likely substantially underestimates the present experimental 
and theoretical uncertainties. Here we also need to input the parton-level polarization, 
P = —0.935. The results of these questionable considerations are shown in Fig. |6| and are 
compared to the CLEO analysis constraints. From this figure we see that the crude estimate 
of the L3 constraints and those obtained above from CLEO are not in conflict and even tend 
to prefer similar regions of the parameter space. The sizes of the two allowed regions are 
rather comparable and substantially overlap. It is also clear from the figure that the L3 data 
certainly excludes both a (V + A) x (V — A) as well as a V x (V — A) interaction as several 
cr as claimed. Before we can draw any stronger conclusions, however, this analysis needs to 
be repeated by L3 themselves with the additional a s and 1/rnl corrections included. We can 
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conclude that future spectra determinations from inclusive decays will indeed be useful in 
probing for RH currents provided high statistics are available and systematic experimental 
uncertainties are under control. Since the L3 analysis is based on a sample of only 10 6 Z's 
it is clear that a higher statistical study can be performed. 
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Figure 6: Comparison of the envelope of the 95% CL allowed regions obtained with CLEO 
data(solid curve) with an estimate of the upper bound obtainable from the analysis of the 
L3 charged lepton and neutrino data(dotted curve). The location of the x 2 minima from the 
L3 analysis is also shown. 
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3 The Left-Right Model and f 



If RH currents do exist, given the CLEO allowed region in the |£| — ca plane shown in 
Fig.l, we want to know if there are any sub-regions of this allowed space that can yield a 
simultaneous lowering of the SM predictions for Bi and n c . To address this question we 
will need to go beyond the physics described by the effective Hamiltonian in Eq.(l) since, 
e.g., we need to discuss non-leptonic decay modes such as b — > cud(s) and b — > ccs(d) and 
the role RH currents may play in these channels. To do this we need to incorporate the 
physics of 7i s i into a larger framework, e.g., the LRM||. We remind the reader that other 
frameworks, such as SUSY loops, compositeness or R— parity violation schemes, are also 
possible sources of effective RH currents. 



In order to be self-contained let us briefly review the relevant parts of the LRM we 
need for our discussion below; for details of the model the reader is referred to ||. The LRM 
is based on the extended gauge group SU(2) L x SU{2) R x £7(1). Due to this extension there 
are both new neutral and charged gauge bosons, Z' , W^, in addition to those present in the 
Standard Model. In this scenario the left- (right-) handed fermions of the SM are assigned 
to doublets under the SU(2)l(r) group and a RH neutrino is introduced. The Higgs fields 
which can directly generate SM fermion masses are thus in bi-doublet representations, i.e., 
they transform as doublets under both SU (2) groups. The LRM is quite robust and possesses 
a large number of free parameters which play an interdependent role in the calculation of 
observables and in the existing constraints on the model resulting from various experiments. 

As far as B physics and the subsequent discussion are concerned there are several 
parameters of direct interest. The most obvious free parameter is the ratio of the SU(2)r 
and SU(2)l gauge couplings 0.55 < k = gn/gi < 2; the lower limit is a model constraint 
while the upper one is simply a naturalness assumption. GUT embedding scenarios generally 
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suggest that k < 1 j31| . For simplicity we will assume that k = 1 in almost all of our discussion 
below. The extended gauge symmetry is broken in two stages. First the SU (2) l x SU (2)r x 
U(l) symmetry is broken down to the SM via the action of Higgs fields that transform either 
as doublets or triplets under SU(2)r. This choice of Higgs representation determines both 
the mass relationship between the Z' and Wr (analogous to the condition that p = 1 in 
the SM with only Higgs doublets and singlets) as well as the nature of neutrino masses; in 
particular, the Higgs triplet choice which we employ here allows for the implementation of 
the see-saw mechanism and yields a heavy RH neutrino. 

After complete symmetry breaking the resulting Wl — Wr mixing is described by two 
parameters, a real mixing angle, 0, and a phase, to. Note that it is usually t = tan0 which 
appears in expressions directly related to observables. The additional phase, as always, can 
be a new source of CP violation. The mixing between Wl and Wr results in the mass 
eigenstates Wip, with a ratio of masses, r = Mf/M%, (with M 2 ~ Mr). In most models 
T is then naturally of order a few times r or less in the large M 2 limit. Of course, W\ is 
the state directly being produced at both the Tevatron and LEPII and is identical to the 
SM W in the <fi — > limit. We note that when is non-zero, W\ no longer couples to a 
purely LH current. Of course if a heavy RH neutrino is indeed realized then the effective 
leptonic current coupling to W\ remains purely LH as far as all low energy experiments 
are concerned. As is well-known, one of the strongest set of 'classical' constraints on this 
model arises from polarized p decay 0, which are trivial to satisfy in the case of a heavy RH 
neutrino and this justifies the appearance of only LH leptonic couplings in Eq.(l). Removal 
of these constraints provides significantly more freedom in the remaining LRM parameter 
space. Thus the tree-level p decay Hamiltonian is just 



W " = ^^Mt (^7AMx)(ex7 A ^ e J , (29) 



A, 
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so that the tree- level definition of Gf is simply 

G F <?I4(1 + rt 2 ) 



y/2 8M 



2 



(30) 



We see that if r and t are of order ~ 10 2 or less the numerical influence of mixing in this 
relationship will be quite small. 

We note that it is important to be reminded that the extended Higgs sector associated 
with both the breaking of the LRM group down to U(l) em and the complete generation of 
fermion masses may also have an important role to play in low energy physics through both 
the existence of complex Yukawa and/or flavor- changing neutral current type couplings. 
However, this sector of the LRM is highly model dependent and is of course quite sensitive 
to the detailed nature of the fermion mass generation problem. For purposes of brevity 
and simplicity and because tree-level neutral Higgs exchange can little influence the decay 
processes we are interested in these too will be ignored in the following discussion and we will 
focus solely on the effects associated with Wig exchange. We do note that these additional 
Higgs fields can potentially play a very important role in loop processes as will be briefly 
discussed later. 

Additional parameters arise in the quark sector; in principle the effective mass ma- 
trices for the SM fermions may be non-hermitian implying that the two matrices involved 
in the bi-unitary transformation needed to diagonalize them will be unrelated. This means 
that the elements of the mixing matrix, Vr, appearing in the RH charged current for quarks 
will be unrelated to the corresponding elements of Vl = Vckm- Vr will then involve 3 new 



angles as well as 6 additional phases all of which are a priori unknown parameters ^] . The 
possibility that Vl and Vr may be unrelated is sometimes overlooked when considering the 
potential impact of the LRM on low energy physics and there has been very little detailed 
exploration of this more general situation due to the rather large parameter space. Certainly 
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as the elements of Vr are allowed to vary the impact of the extended gauge sector on B 
physics in general will be greatly effected. 



Some well-known constraints on the LRM, such as Tevatron direct W searches 



are quite sensitive to variations in VrP^I as wen as the properties of the RH neutrino and 
W 2 masses as low as 450 — 500 GeV can very easily be accommodated by the present data. 
To be conservative, and with future Tevatron searches in mind, however, we will assume 



below that M% > 720 GeV|33[], i.e., r < 0.012, for any Vr implying that the magnitude 



of t is also less than ~ 0.012. Other constraints on the LRM parameter space involve 



loop processes such as the Ki — K$ mass difference fl35|, |3ql and b — > 57^3. Clearly the 
bounds obtained from these processes depend not only on the gauge sector but also on all 
the particles that can participate in the loops such as SUSY partners, extra Higgs fields, 
additional heavy fermions, etc., whose existence is sensitive to the finer details of the model. 
These possibilities are beyond the necessities of the current discussion where we are solely 
interested in tree-level B decays. Our philosophy as outlined in the introduction will be to 
leave for now all discussions of loop graphs which display any sensitivity to the details of the 
LRM spectrum and take these issues up briefly later. 

Using the definitions above for the LRM parameters we can now express £ in terms 
of these more fundamental quantities; we find 

(31) 



;i + re-) v e 



Note that we can absorb the sign of t into the phase u here so that t can be treated as 
a positive parameter in our discussion below. As mentioned above we will take k = 1 
for simplicity in our numerical analysis below; to first order it simply rescales the value 
of t. Employing the results from BurasfllH] for the value of Vj£ mc (_D*) from inclusive B 
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semileptonic decays, we can invert the expression above to obtain 



|^H(39.9 ± 2.2). 10 -3 [1 + |e ^ £01, (32) 



so that for typical values such as |£| = 0.2, ca=0, and x = 0.29, we obtain 

10- 



V£\ = (0.782 ±0.042) 



Kt 



[l + 0(r,rt 2 )], (33) 



which suggests that \V£\ is reasonably large and perhaps of order unity over most of the 
allowed parameter space shown in Fig. |]. From these considerations we learn several things 
which follow immediately from the unitarity of Vr: (i) A large value for \V^\ implies that 
the sum | V^\ 2 + | V c f | 2 is small thus somewhat suppressing potential RH contributions to the 
decays b — > ccs(d), which is fortunate for charm counting purposes. If either of these elements 
were large one might expect a significant increase in n c due to RH current contributions. As 
we will see below, just the opposite occurs. As will be noted, this also assists in suppressing 
RH contributions to K L — K s mixing, (ii) Since unitarity requires \V^ d \ 2 + \V^\ 2 + \V^\ 2 = 
\ v ^\ 2 + \ v cb I 2 + \ v i I 2 ^ follows immediately that \V^ d \ 2 + \ V^ S \ 2 > \V C %\ 2 . However, since 
\V c b | 2 is apparently large this inequality implies that the sum \Vj^\ 2 + \ Vj^\ 2 is larger still. 
This would mean that decay modes such as b —>■ cud(s) may receive large RH contributions. 
We note that if we further assume that \Vj^\ 2 « \Vj^\ 2 these RH contributions may also 



lead to an increase in K production in B decays |J8fl, which it has been argued is a signal for 



enhanced b — > sg. Also if \V^\ 2 « \V,^ S \ 2 one finds that the Tevatron search reach [ 33 1 for 



W 2 would be seriously degraded by about a factor of 20] in mass. (Hi) It would appear that 
\VJI\ will be too small to significantly influence b — > u processes, though this needs further 
examination, (iv) A large implies that the sum \V^\ 2 + \V t f \ 2 is also large PD|] with 
implications for the complete structure of Vr that we will ignore for now but will return 
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to haunt us in our discussion below, (v) The fact that unitarity requires \V£\ < 1 itself 
provides an additional constraint on the remaining LRM parameters. 



4 Non-leptonic b —> c Decays with RH Currents 



As a final step in our analysis we need the complete non-leptonic Hamiltonian; at the tree- 
level this can now be written down immediately. For the sample case of b — > cud we can 
write, following the notation in Refs-ji^j |39|, 
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(34) 



where O2L = (c7 a1 Pl&)(^7 a '-Pl^), Oul = (c r y fJi PRb)(d'y fJ, PLu), etc. and where Pl,r are helicity 
projection operators. At the weak scale the operator coefficients are given by 
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(35) 



Note that if we neglect the light quark masses the appropriate phase space functions for this 
particular decay mode will be given by / and g. The modifications necessary for the study 
of the decay b —* cus are obvious. Similarly for the corresponding decays b — > ccs(d) we 
simply change the appropriate CKM factors in the above and employ the appropriate phase 
space functions, f c and g c , which are given by the phase space integrals I± 2 in Section 2 with 
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the replacement y — > x. For x = 0.29 these are found numerically to be f c ~ 0.222 and 
g c ~ —0.086. The neglect of the strange quark mass, m s ~ 100 — 150 MeV, is found to be 
an excellent approximation here. 

To proceed with this calculation we need to compute the QCD corrections associated 
with the Renormalization Group running from the weak matching scale down to /i ~ m^. 



To this end we follow the analysis of Bagan et a/.|4(J which allows us to write the partial 
width for this process as 



T(b -> cud(s)) = T S M 



1 + rji + 7/2 + r/ 3 



(36) 



where T SM = 3XiT f\V cb \ 2 (\V^ d \ 2 + \V^ S \ 2 ), T is the canonical \x decay width with the 
replacement fi — > m^, and X\ represents the results of SM QCD corrections(to which we will 
return below). The rji are LRM contributions which given by 
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t(l-r) 



K| 2 + 



\d\ 2 y- 
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1 + 



K 2 (r + t 2 ) 
(1 + rt 2 ) 



withy = {\V u R d \ 2 + \V u R s \ 2 )/{\V u L d \ 2 + \V u L s \ 2 ) ~ (|K5| 2 +|K2| 2 ). As pointed out in the discussion 
above, if |V^| is large we anticipate that y is near unity. For the decays b — > ccs(o?) we make 
the obvious CKM replacements and the change the Xj — > X- , f,g—> f c , 9c and let y — >• y c 
where y c = ( | V£ \ 2 + 1 V£ \ 2 ) /( | | 2 + 1 V£ \ 2 ) ~ 1 - | V£ \ 2 , with the last near equality resulting 
from unitarity and the fact that \V^ b \ 2 is very small. If \V C R \ 2 is large then clearly y c must 
then be small. 
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At leading order(LO) in QCD the Xi = X[ are completely calculable and are simple 
polynomials in the parameter 



a s (M } 



w , 



3/23 



(38) 



and its inverse; here we will assume that a s {Mz) — 0.118 and fi ~ m^. Explicitly, we obtain 



X, 



X, 



2z 4 + z' 



\z 2 + z~ m 



(39) 
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Az A + Az~ A + 2z~ b - z 



-12 



where we have made use of the results of Altarelli and Maiani|13| as well as Cho and 
Misiak|57|]. Note all Xi — > 1 as z — > 1 and the QCD corrections vanish. In the SM, NLO 
multiplicative corrections to the LO values of X\ and X[ are now known p^] to be ~ 1.061 
and ~ 1.29, respectively, for /i = m b , x = 0.29, and using pole quark masses(m b = 4.8 GeV), 
both of which we adopt in the numerical analysis below. Unfortunately, the correspond- 
ing NLO corrections to X 2 ^ and X' 2 3 are not yet known. The best we can do until such 
calculations are performed is to follow Voloshin's philosophy and assume the multiplicative 
corrections in these cases are essentially the same as those for X\ and X[. Since, as we 
will see below, we will be more interested in the shifts in the values of n c and Bi due to 
RH currents than the values themselves, we anticipate that this assumption may be a fair 
approximation. We note that in making this assumption we are also ignoring the possibility 
that the detailed LRM particle spectrum may lead to substantial modifications in these SM 
values, in particular, those contributions arising from penguins. These assumptions need to 
be verified by future direct calculations. 
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5 5B^ and 5n c 



From the discussion in the previous section we are ready to calculate both 5 Be = Bg(LRM) — 
Bg(SM) and 5n c = n c (LRM) — n c (SM) where the SM results are given by the above 
expressions in the limit where all RH couplings are turned off. As is well known, the combined 
experimental and theoretical situation is quite puzzling. From the reviews of both Drell and 
Sachrajda0 we see that B t = 0.1018 ± 0.0040 on the T(4S) while B t = 0.1095 ± 0.0032 at 
the Z. Similarly, n c = 1.119 ± 0.053 and 1.202 ± 0.067 on the T(4S) and Z, respectively. 
Numerically, in the SM limit our calculations essentially reproduce the earlier results of 
Bagan et at. [fHi ] which we have closely followed; in this limit we obtain Be = 0.123 and 



n c = 1.24 for the SM predictions assuming x = 0.29 and fi = nib- We will implicitly assume 
that there are no new b —>■ no charm final states, such as b — > sg, which are enhanced due to 
RH currents. It is clear that if we take these experimental numbers at face value we would 
like to decrease the theoretical predictions for Be by 0.015 — 0.020 and n c by at least 0.03. 

Our analysis consists of an extensive scan of the model parameter space spanned by 
r, t, |£|, ca and y and demanding that a number of requirements be satisfied simultaneously. 
Our input parameters are chosen as follows. We begin by picking a 'point' inside of the CLEO 
allowed region in the ca — |£| plane so that this constraint is already satisfied. We assume 
the scale size of the ca — |£| grid to be 0.01 x 0.01 so that are approximately 1.5 • 10 4 points 
in this sample. Next, we choose a value for the two LRM parameters r and t; for simplicity 
k is set to unity. Keeping in mind the CDF/DO bounds and the strong suggestion that 
t cannot be much larger than r, we let r = 0.0025, 0.005, 0.0075, 0.010 or 0.012 and allow 
t to vary over the range to 0.012 in steps of 0.0005. (Remember that due to the phase 
freedom in angle A we can treat t > in this discussion.) Clearly if the Wi mass is too large 
and/or the mixing angle is too small the effects of RH currents will not be of a noticeable 
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magnitude. This restricts our attention to W 2 masses in the approximate range 730 — 1600 
GeV. Thus we see that for every choice of ca and |£| there are 120 pairs of (r, t) values giving 
us a total of ~ 1.8 • 10 6 points to examine in the r — t — |£| — c A parameter subspace. 

The first constraint we impose is the requirement that \V£\ be less than unity by 
using Eq. (32). Of course if this constraint is not satisfied for any of the r or t values this 
point on the ca — |£| grid is removed from any further consideration. If satisfied, the result 
fixes the value of y c in the subsequent calculations. To proceed we must choose a value of 
y in the range < y < 1 which we do in grid steps of 0.01. We then impose our second 
constraint that y > \V^\ 2 so that only the larger of the y values survive. Out of the original 
~ 1.8 ■ 10 8 points in the five-dimensional r — t — |£| — ca — y parameter space being scanned, 
only ~ 27.5 • 10 6 survive these first two constraints. 

For these remaining points we next calculate 5Bp and 5n c for each particular choice 
of input parameters and impose our final loose requirement that 5Bp < —0.01 and 5n c < 
—0.025. Again, if these constraints cannot be met at a particular point on the ca — |£| 
grid, independently of the chosen values of r, t and y, it is removed. Only 6284 points in 
the the r — t — \^\ — — y five-dimensional parameter space now remain; this number is 
further reduced to 972 if we strengthen our requirement on 5n c to be < —0.03. It is clear 
from these numbers that a rather high degree of fine-tuning is required to push Bp and n c 
in the proper direction and to produce shifts with interesting magnitudes. For most values 
of the parameters the resulting shifts in Bp and/or n c are much too small to be of interest. 
The combination of these requirements is found to be extremely demanding on the model 
parameter space, yet two distinct sub-regions do survive. 

If one plots the values of 5Bg and 5n c for the survivors we find that they essentially 
lie only along two straight lines in the 5Be — Sn c plane with the choice of line depending 
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Figure 7: Location of surviving points in the 8B^ — 5n c plane. The 972 survivors of the 
Sn c < —0.03 cut are shown explicitly. The lines represent smoothed versions of the actual 
locations. The solid(dash-dotted) line corresponds to the solutions with ca > (<)0. 
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upon the sign of ca as shown in Fig. [?|. The corresponding location of these same points 
with 5n c < —0.03 projected onto the ca — |£| plane are shown in Fig. [8]. It is amusing to note 
that the points with ca > lie within the region associated with the fit to CLEO's data on 
the x distribution in B — ► D*{-^ Dii)lv obtained above. Assuming 5n c < — 0.025(— 0.03) 
approximately 92.5(75.1)% of the survivors are found to lie in the ca > region. The 
fractional volume of the 5n c < —0.025 parameter space which also allows 5n c < —0.03 is 
~ 15.5%. While the ca < parameter space is only reduced to 51.4% of its previous size by 
strengthening this Sn c cut, the ca > subspace is drastically reduced to only 12.6% of its 
previous population by this same cut. 

What are some of the various properties of the parameter space points that satisfy 
all our requirements? Mostly they are exactly what one would naively expect. First, all of 
the 972 survivors have t > 0.0095 since larger mixing angles are required to enhance the 
contributions of the RH currents. Second, in all cases \V£\ > 0.908 and there is a significant 
preference for larger values of r, i.e., there are only 4(37) cases with r = 0.0025(0.005). 

6 Discussion and Conclusions 

The chirality of the b — > c coupling is one of the most important quantities in B physics. 
The original work of Gronau and Wakaizumi demonstrated to us just how little was actually 
known about this coupling at that time. Since then, after extensive theoretical and exper- 
imental effort, the situation remains far from being completely clarified. While CLEO and 
ALEPH have certainly demonstrated that the h — > c coupling is dominantly LH in agreement 
with the SM, their results remain consistent with the possibility of a sizeable RH coupling. 
Furthermore, the interpretation of the low value of the A& polarization obtained by L3 re- 
mains ambiguous and could either be a first signal for RH currents or simply a sign of our 
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Figure 8: Locations of the zones containing the 972 surviving points in the ca — |£| plane, 
in comparison to the envelope of that allowed by CLEO at 95% CL, which simultaneously 
satisfy Sn c < -0.03 and SB e < -0.01. 
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ignorance of the strong interactions. All of these experimental analyses have been based on 
relatively small sample sizes and need to be repeated and improved upon. 

As we saw in the analysis above, the exclusive B — > D* semileptonic decay provides us 
with a large number of observables that can be used to probe for RH couplings of reasonably 
small strength. In addition to the overall partial width, expressible in terms of exc , the 
measurements of the cos dp and cos^y distributions lead directly to the quantities Afb and 
Tl/Tt, respectively. By using HQET we performed a fit to the present CLEO results for 
these quantities and demonstrated that the current bound on the RH coupling strength still 
remains rather poor especially if it is allowed to be complex. Improved statistics available 
at upcoming B factories will help tremendously here. Furthermore, while we showed that 
the q 2 distribution was not very sensitive to RH couplings, the \ distribution was found 
to be particularly so and yielded tantalizing indications for the existence of RH currents. 
Present CLEO data was shown to indicate that future measurements of this distribution will 
be extremely useful in either constraining or discovering RH couplings. More recent but yet 

unpublished CLEO results [f20j seem to strengthen the case for the existence of right-handed 
currents based on the x distribution. 

The low Afe polarization result obtained by ALEPH also remains tantalizing and 
certainly needs updating. Unfortunately, given our incomplete knowledge of QCD, any in- 
terpretation of the result in terms of RH currents can not be made at present. However, 
if high precision measurements of the lepton and missing energy spectra become available 
with only a factor of a few increase in statistics, we saw in the analysis above that suffi- 
cient observables do exist to separate the two possible explanations. Further observables 
may be found to strengthen any conclusions one may draw from future data. The low A5 
measurement seems to be confirmed by as yet unpublished data from both ALEPH and 
DELPHI [0. 
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Under the assumption that b — > c RH currents do exist consistent with the bounds 
from CLEO, we have tried to address the question raised by Voloshin as to whether such 
new interactions could assist in solving the long standing problem associated with Bp and 
n c . To address this point we needed to go beyond the model independent results of the 
previous section and incorporate our b — > c RH coupling scenario into a larger framework, 
the most natural one being the LRM. Within this scheme, making a number of assumptions 
about both the detailed particle spectrum of the model and the nature of the NLO QCD 
corrections to the RH pieces of the nonleptonic Hamiltonian operator coefficients, we were 
able to demonstrate that two small regions of the LRM parameter space do exist that push 
both Bp and n c in the right directions with sufficient magnitudes to be phenomenologically 
interesting. These small parameter space regions result from a reasonably highly tuned set 
of LRM parameters and in all cases Vjjf was found to have a magnitude of order unity. 

Hopefully measurements at the new B factories which are soon to turn on will yield 
signals for physics beyond the Standard Model. Perhaps right-handed currents will be among 
them. 
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APPENDIX 



In this Appendix, we outline some of the implications of the scenarios discussed above 
that lead to lower values of both Bi and n — c while satisfying the CLEO constraints. Such 
results, for example, may lead one to speculate on just what form the matrix Vr might take 
if this type of solution to the Bg_ — n c problem were to be realized. This will directly lead to 
a number of wide ranging implications in all low energy sectors of the theory and not just 
in B physics. (In fact, there are too many for us to comment upon here with any depth 
of discussion.) Unfortunately, we do not yet have available a global analysis of RH current 
phenomenology for arbitrary forms of Vr with a left-right mixing at the per cent level. Such 
an analysis would be extremely useful for our discussion but is far beyond the scope of the 
present paper. 



If we hypothesize [j36[ [39f that in each row or column there is a single element with a 
magnitude near unity, as is true for the conventional CKM matrix, then there are only two 
RH mixing matrices which allow for large V£. Following the notation employed in our earlier 
work[35], we can write these 'large element' forms symbolically, neglecting any phases, as 



^ 1 ^ 



M 



c 



1 
1 



^ 1 ^ 



M 



D 



1 

1 



(40) 



with the true Vr being a perturbation about one of these skeletons, just as the CKM is 
a perturbation about the diagonal unit matrix. As noted elsewhere [j39j] the structure of 
these matrices combined with small values of t allow us to easily circumvent the traditional 
constraints on the LRM from the magnitude of Kl — K$ mixing. However, we also observe 
here the necessity that at least one of V t ^ or V t f is large which can lead us to some potential 
problems with the observed rate for and limits upon the processes b — > d, S7 ||30|| . As has been 
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discussed in previous analyses of the b — > 57 process within the context of the LRMpT], [L2 



|39| , interference terms between the RH and LH W\ — W2 contributions obtain enhancements 
by helicity-flip factors of mt/rrib though they are also simultaneously suppressed by a factor 
of t. While the pure SM piece is proportional to the product V t ^V t ^ s , these new interference 
terms are correspondingly proportional to either V t ^V t ^ s or V t ^V t ^ s ; the later being quite 
small in our case. However, here we have seen that at least one of the products V t ^V t ^ s is 
of order unity This implies that in the models we have found here such LH-RH interference 
terms arising from W12 ~ t quark penguins may be dangerously large, by factors of order 
10, in at least one of the b — > 57 or b — > modes. Of course one may argue that we are 
most likely quite ignorant of the true LRM model spectrum and that loop contributions 
from the non-Wj — t diagrams may eliminate this problem. As is well known, SUSY and 
charged Higgs exchanges can, for example, yield significant contributions to these penguins 
and can possibly leading to a fine-tuned cancellation amongst the various pieces. This is an 
unnatural, yet potentially possible, solution. 

Another potentially important constraint arises from the determination of the relative 
branching fractions for the B — > ipir and B ipK modes by CLEO[^l| to be 4.3 ± 2.3%. 



This measurement was inconsistent with the predictions of the Gronau and Wakaizumi model 
which predicted a ratio of ~ 10~ 7 . For the class of models presented here, this result implies 
only that Vj$ is probably somewhat smaller than V^, which is not an unexpected result. 

As is well-known there are many other non— B physics constraints on the form of Vr 
which need to be examined. These are mostly concerned with the specific elements Vj^ s ; 
some of these have a rather long, even controversial, history. Many of these constraints have 
been extensively reviewed in detail some time ago by Langacker and Uma Sankar|RB| and 



as stated above it is beyond the scope of the present paper to discuss them at any length 
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except for several comments. These low-energy constraints include, amongst others, potential 



violations of CKM universality, as suggested by Wolfenstein||42||, and/or violations of PCAC 



relations, as suggested by Donoghue and Holstein|EJ. For the universality constraint, we 

note that Buras@ reports Ei \V% e// | 2 = 0.9972 ± 0.0013, which is more than 2a below 
the SM expectation, perhaps hinting at new physics, //no other new physics sources enter 
other than the existence of Vr, we can use the result of this sum to constrain both -Re(V^J 

even when nt ~ 0.01. For example, using \V^ d \ ~ 0.98 and \V^ S /V^ d \ ~ 0.22, this constraint 
implies 



0.133 ±0.066 



10 



-2 



Krf|cosA d + 0.221 V^l cos A s 



(41) 



where Aj is sum of u plus the phase of V^. This constraint is easily satisfied for either of 
the two forms of Vr suggested above assuming reasonable Aj. Some possibilities, suggested 
by an earlier analysis of Matrix D[]39| are to either have V^ s essentially unit magnitude but 
with a rather large phase together with \V^\ ~ A 2 ~ 0.05 with arbitrary phase or to have 



instead \V^ d \ ~ A ~ 0.2 with bot Vj^ s having sizeable phases. (At this point we remind 
the reader, however, that in extended gauge theories such as the LRM there can be other 
potentially significant contributions to universality violation, e.g., Z' exchange, as has been 



discussed by Marciano and Sirlin|44j].) Interestingly, such possible solutions are also found 
to easily satisfy a number of additional constraints including those from PCAC{|3| (though 
these need to be updated), those from muon capture on 3 He[[15|], and those on the phase of 
9 a/ ' 9v i n neutron beta decay H9]. Similarly, the scaling of the strengths of the V, A currents 



imply that the extracted value of the ratio [sa) 9v)l f-K is exactly the same as in the SM 
with no violation of the Goldberger-Treiman relation^] occurring in the presence of RH 
currents. While safely avoiding all these bounds, however, these solutions do not help in 
explaining a possible disparity between the values of V^ d extracted from neutron decay and 
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that obtained from + — > + and 19 Ne beta decay RIB |. This at the very least would require 
a quite sizeable V^ d . 

In these same scenarios one might expect somewhat larger effects due to RH currents 
to now appear in the strange quark sector. Perhaps one of the most significant effects of 
RH currents here, apart from overall changes in normalizations of constants, is in the F and 
D parameters describing hyperon decay. The values extracted for these parameters from 
AS = and AS* = 1 transitions, corrected for the not yet completely understood £{7(3) 
breaking effects, would appear somewhat different. The reason here is clear: the ratio of 
the axial-vector to vector coupling constants in the two cases are shifted away from their 
SM values by different amounts depending on the form of Vr. Although the data remains 
rather poor, this possibility is not unsupported by the recent analysis of Ratcliffe pH . The 
implications are, of course, far reaching and extend as far as tests of the Bjorken Sum 
Rule[|50|. We also remind the reader of the well known |fl6f| potential discrepancy between 



the value of V us extracted from the vector current coupling in K e3 decays and that from 
hyperon decay data, which probes both axial-vector as well as vector couplings. 

Clearly, is a possible signature of RH currents arises in B decays, the search for their 
influence elsewhere becomes ever more important. 
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